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Results of simulations of flow of an arc-heated stream around a 14-inch diameter 45° 
sphere-cone configuration are presented. Computations are first benchmarked against 
pressure and heat flux measurements made using copper slug calorimeters of different 
shapes and sizes. The influence of catalycity of copper on computed results is investigated. 
Good agreements between predictions and measurements are obtained by assuming the 
copper slug to be partially catalytic to atomic recombination. With total enthalpy estimates 
obtained from these preliminary computations, calculations are then performed for the test 
article, with the nozzle and test article considered as an integrated whole - the same 
procedure adopted for calorimeter simulations. The resulting heat fluxes at select points on 
the test article (points at which fully instrumented plugs were placed) are used in material 
thermal response code calculations. Predicted time histories of temperature are compared 
against thermocouple data from the instrumented plugs, and recession determined. Good 
agreement is obtained for in-depth thermocouples. 
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Subscripts 
Air = 

Air+ 

Ar 

arc = 

bulk 

cal = 

CL 

e = 
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Acronyms 

AHF 

ARC 

CFD 

DS2 

DAS 

EB2 

HEAT 

IHF 

LSAT 

MSL 

OME 

PICA 

RTV 

SPORE = 
SPRITE = 
TTT 


primary flow stream 

secondary flow stream (“add air”) in the plenum 

argon 

are eolumn 

bulk quantity in the plenum (i.e., mass-averaged value) 

ealorimeter 

eenterline 

boundary-layer edge 
hemisphere 
stagnation point 
nozzle throat 
wall 

aerodynamie heating faeility 

Ames Researeh Center 

eomputational fluid dynamies 

Deep Spaee 2 

data aequisition system 

enthalpy by energy balanee 

Hollow aEro thermal Ablation and Temperature 

interaetion heating faeility 

large-seale artiele test 

Mars Seienee Laboratory 

outer mold line 

phenolie impregnated earbon ablator 

room temperature vuleanizing 

Small Probes for Orbital Return of Experiments 

Small Probe Re-entry Investigation of TPS Engineering 

through the thiekness 


I. Introduction 

I n an effort to fulfill a ‘test what you fly’ paradigm/’^ a 14-ineh diameter entry probe named ^Small Probe Pe-entry 
investigation of TPS Engineering (SPRITE) has been developed at NASA Ames Researeh Center. It is envisaged 
that this elass of small probes, with suitably tailored trajeetories, will serve as robotie flight test beds for TPS 
materials^ or perhaps even other vehiele subsystems (the SPRITE design is quite flexible). With some modifieations, 
the same small probe design eould also be used to return biologieal samples from orbit, whieh is the prime objeetive 
of the Pmall Probes for Orbital Return of Experiments (SPORE) projeet led out of Georgia Institute of Teehnology."^ 
As a first step towards establishing the feasibility of the ‘test what you fly’ paradigm, a full-seale (14-ineh 
diameter) model of SPRITE was tested in the I8-ineh nozzle of the Aerodynamie Heating Faeility (AHF), an are jet, 
at NASA Ames Researeh Center. This series of are jet tests for SPRITE was a follow-on to the LSAT or Large- 
Seale Artiele Test^ performed in the NASA ARC 60 MW Interaetion Heating Faeility (IHF) in late 2008/early 2009, 
and the test of a full-seale red oak model of DS2 or Deep Spaee-2^ in the NASA ARC 20 MW Aerodynamie 
Heating Faeility (AHF) in late 2009/early 2010. There were no issues with mass eapture by the diffusers in either 
faeility (AHF or IHF) for blunt bodies up to 15 inehes in diameter. Of eourse, large nozzles (18- and 21-ineh exit 
diameter) had to be used in testing artieles of this size. 

Building on this initial sueeess with large-seale artieles, two identieal models (SPRITE-Tl-1 and SPRITE-Tl-2 - 
T1 indieating the ehoiee of baek shell geometry) were fabrieated and tested in the AHF with the 18-ineh nozzle 
installed. Both these test artieles were 14 inehes in diameter and had a 45° sphere-eone forebody (like DS2) made 
of PICA (Phenolie Impregnated Carbon Ablator) bonded on to a 1/8^^-ineh thiek aluminum shell using RTV 
adhesive. The aft portion of the test artiele was a eonieal frustum (15° eone angle) with LI-2200 bonded on to the 
aluminum shell. The PICA shell had a thiekness of 1 ineh over the frustum and 1.36-ineh thiek on the spherieal nose 
eap, while the LI-2200 tiles on the aft shell were 0.7-ineh thiek. Eaeh model was fully instrumented with (a) 3 in- 
depth thermoeouples and a HEAT sensor^ imbedded in 1.3 -ineh diameter plugs in the heat shield (similar to those 
used on the heatshield of the Mars Seienee Laboratory or MSL^), (b) several thermoeouples bonded to the aluminum 
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substructure, and (c) a few strain gauges. Data from the strain gauges and some of the thermoeouples and gauges 
were aequired by the on-board data aequisition system (DAS),^^ while data from the other thermoeouples were 
routed to the faeility-provided DAS, thereby enabling a eross eheek on the in situ measurement eapability. Apart 
from the primary objeetive of testing flight-seale artieles in an are jet, the experiments also served the purpose of 
demonstrating the feasibility of in situ measurements of temperature, strain and reeession using a data aequisition 
system integrated internally with the test artiele. 

Modeling and simulation have been an integral part of the SPRITE are jet experiments, and were used in test 
planning and post-test analysis. For instanee, pre-test flow field eomputations were performed to assess flow 
bloekage issues (if any).^ Results from these eomputations were used in estimating the time of exposure of the test 
artiele to the are-heated stream so that the internal temperature of the probe would not exeeed the eonstraint of 
battery operation of the internal data aequisition system. The eombination of simulations tools - primarily Dplr^^, 
Fiat^^, and Mar&^ - were then used in post-test analyses in order to verify/validate these tools against aetual 
measurements. 

The present paper also provides details of the proeess used in simulation of the flow of an are-heated stream 
expanded through a eonieal eonvergent-divergent nozzle over a test body (ealorimeter or the test artiele). The 
primary objeetives of the DPZi?-based simulation proeess^^ are: (1) to assess if there are bloekage issues with 
eandidate test geometries, (2) to provide initial heating distributions over the test artiele so that thermal -struetural 
analyses ean be performed to estimate the required exposure time, (3) to provide post- test heating estimates, based 
on as-run eonditions of the are heater, for final thermal- struetural analysis, and (4) to determine the thermal response 
(in-depth thermal environments and/or reeession, if any) of the materials used on the fore and aft shells of the test 
geometry. In addition, these are jet tests eould faeilitate a better understanding of the surfaee eatalyeity of eopper, 
the material used in ealorimetry. 


II. Test Geometry 

The test geometry is a 45° sphere-eone of 14 inehes base diameter {i.e., 7?/, = 14 in). The geometry is shown in a 
eoordinate system eentered on the nozzle exit plane (NEP) (Fig. 1). The nose radius, is 4 inehes, and the radii of 
the two shoulders, Rs and Rc (fore and aft), are 0.4 and 0.25 ineh, respeetively. The eonieal frustum of the aft shell 
has a 15° inelination to the horizontal, and the total axial length, of the test artiele is 11 inehes. Also shown in 
Fig. 1 are the loeations of the instrumented plugs used in the test. The loeations of the plugs are approximate and 
measured from the tangeney points of the front eonieal frustum and the shoulder torus. For the purpose of CFD 
analysis, i.e., predietion of surfaee aerothermal environments (pressure, heat flux, and shear), only the outer mold 
line (OMF) of the test artiele is neeessary; the time-varying distribution of heat through the thiekness of the 
materials that make up the test artiele is determined by a materials thermal response eode sueh as FlAT^^ or a 
thermal-struetural analysis eode sueh as Marc}^ with the CFD predietion of environments as boundary eonditions. 



Figure 1. SPRITE-Tl geometry tested in the 18-inch nozzle of the 20 MW Aerodynamic Heating Facility 
(AHF). The apex of the model is located 12 inches from the nozzle exit plane (NEP). 
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III. Pre-Test Computations 

Testing at large seale (12 to 14-ineli diameter test artieles) provides meaningful spatial development of a 
boundary layer on the wetted aereage as eompared to the 4-ineh diameter test artieles (usually simple shapes like 
flat-faeed eylinders, or hemispheres, or iso-q) that form the basis of typieal stagnation point testing in the ARC are- 
heated faeilities. However, the faeility dimensions (nozzle and diffuser) plaee an upper bound on the projeeted 
diameter of the test artiele, with flow bloekage, i.e., inability of the diffuser to eapture the mass flow, being the 
primary eonsideration. 

A simple way of establishing the model size limit is to build wooden models (of different sizes) and expose them 
to an are-heated stream. Sueh a proeedure poses some risk to the faeility. A simpler way than the build-and-test 
approaeh, and a quieker one, is to run a number of CFD simulations with models of different shapes/sizes, and 
thereby determine the maximum size for whieh the diffuser is still able to eapture the flow. A single wooden repliea 
of the final shape then needs to be exposed to the are-heated stream to establish that testing at that seale ean be 
aeeomplished safely in the faeility. This simpler CFD-based approaeh was ehosen in the design of the sueeessful 
LSAT^ eondueted in the 21.5-nozzle of the IHF, and the same approaeh was taken for the design of the SPRITE test 
artiele as well. The SPRITE test artiele was designed for the 18-ineh nozzle of the AHF. 

The CFD proeess, built around v3.05 of Dplr,^^ is doeumented in the paper of Prabhu et al}^ Briefly, the 
simulation proeess treats the problem as a gas dynamie one and exeludes the eonstrietor are heater. The are-heated 
gas mixture of air and argon is assumed to be in thermoehemieal equilibrium at the inlet to the eonvergent-divergent 
nozzle, and inflow profiles of speeies mass densities, temperature, and veloeities are estimated using an in-house 
software utility ealled Nozzle _Thro AT JCONDITIONS,^^ whieh uses Cea^^ as the equilibrium eomputational engine. 
This utility requires the flow rates of air and argon, the bulk enthalpy of the are-heated gas mixture, and an estimate 
of the Maeh number at the inlet. Computations are then performed with the inflow profiles preseribed as pointwise 
boundary eonditions in Dplr. The eomputational domain ineludes the eonvergent-divergent nozzle, the test artiele, 
and part of the diffuser. Sinee the methodology eonsiders flow non-uniformities (if any) to oeeur in the radial 
direetion only, the axisymmetrie version of the CFD proeess suffiees for axisymmetrie test artieles. 

In the pre-test eomputations, the SPRITE geometry was identieal to the DS2 flight vehiele - the 45° sphere-eone 
forebody had a base diameter of 13.8 inehes, and the aft shell was a seetion of a sphere (required for stability in the 
low supersonie through subsonie portion of flight). Caleulations were performed for the AHF operating at eonditions 
elose to its maximum - a eurrent of 2000 A and an air-argon total mass flow rate of 360 gm/s. At these eonditions, 
the are eolumn develops a pressure of 814 kPa. Figure 2 shows piteh plane eontours of Maeh number obtained from 
the eomputations. Clearly, the diffuser eaptures the test artiele bow shoek, and there are no issues related to flow 
bloekage for a 1 4-ineh diameter test artiele. It should be noted here that bloekage is not an issue beeause the 45° 
sphere-eone forebody sweeps the shoek baek. For sphere-eone geometries with larger eone angle, there eould be a 
potential issue with a 1 4-ineh diameter body. 
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Figure 2. Pitch plane contours of Mach number for the 18-inch nozzle of the AHF with a 14-inch diameter 45° 
sphere-cone geometry (DS2) placed at 12 inches from the nozzle exit plane. 

Since computations indicated that there are no problems with testing a 14-inch diameter geometry in the 18-inch 
nozzle of the AHF, a red oak model of the same geometry was constructed and exposed to an arc-heated stream. The 
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current and mass flow rates were identical to those used in the computations. The left half of Fig. 3 shows a still 
photographic image of the wooden artiele in the arc heated stream, and the right half of the same figure shows the 
eontours of CFD-predicted statie pressure around the test artiele along with stream traees. 



Figure 3. A still photographic image of a red oak replica of DS2 (13.7-inch diameter) in the arc-heated stream 
of the 18-inch nozzle of the AHF. The image on the right has pitch plane contours of predicted static pressure 
along with stream traces. 


The dark areas in the still photograph are the regions where the wood is eharred by the are-heated stream. It is 
interesting to note that some portion of the aft shell is not charred at all. The reason for this ean be gleaned from the 
image on the right half of Fig. 3. The flow over the aft shell remains attaehed for some distanee and eventually 
separates. The region where the wall-bounded shear layer remains attaehed to the aft shell is the one that is eharred. 
Post-test examination of the test artiele showed good eorrespondence between the length of the charred region and 
length of the attaehed flow predieted by CFD, giving some confidence in the predietive capability, even though 
quantitative information is not available in the wooden model test. Furthermore, the “similarities” between the 
luminous regions of the still photograph, and features in the predietion (shear layer angle, nozzle shear layer-shock 
interaction) add to this confidence. 



Axial coordinate (from NEP), in 



10 15 20 25 30 

Axial coordinate (from NEP), in 


(a) Pitch plane - entire computational domain (b) Test article 

Figure 4. Pitch plane contours of Mach number for the 18-inch nozzle of the AHF with a 14-inch diameter 45° 
sphere-cone geometry (SPRITE-Tl) placed at 12 inches from the nozzle exit plane. Also shown is a magnified 
view of the test article with contours of static pressure and stream traces in the pitch plane. 


In order to avoid difficulties in machining LI-2200 (Shuttle tiles) for an aft shell with eurvature, as in the DS2 
shape, some simplifieations were made for SPRITE-Tl. The rounded aft shell was replaced by a eonical frustum 
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(cone angle of 15° to the horizontal), and CFD computations were performed again for the new geometry. The flow 
eonditions were the same as those used in the DS2 ealeulations. The results of the new computations are shown in 
Fig. 4. Pitch plane contours of Mach number in the entire computational domain are shown in Fig. 4a. Clearly, as in 
the DS2 case, SPRITE-Tl does not have any flow blockage issues. A magnified view of the test artiele with pitch 
plane eontours of statie pressure and stream traces is shown in Fig. 4b. The forebody flow for SPRITE-Tl is almost 
identieal to that of the DS2 (Fig. 3b). The flow on the aft shell remains attaehed over the entire aereage, separating at 
the aft shoulder. The predicted surfaee pressure, heat flux, and shear distributions were used as initial conditions for 
detailed thermal- structural computations. These thermal- structural computations were performed to determine the 
rate of heat soak into the structure and determine a priori an estimate of the time of exposure for the actual arc jet 
model before the internal temperature exceeded a pre-set value. More details of the use of CFD predictions in 
thermal-structural computations can be found in the paper of Agrawal et al}^ 

IV. Computational Analysis and Results 

Two fully instrumented models of SPRITE-Tl were tested in the AHF (Entry 295) at NASA ARC. Flow field 
eomputations were neeessary to reconstruct the pressure and heat flux environments experienced by the SPRITE 
article at actual test eonditions. The axisymmetric simulation methodology developed by Prabhu et al}^ for the arc 
jets at NASA ARC requires as inputs (a) the mass flow rates of air and argon, and (b) the specifie total enthalpy on 
the nozzle eenterline. The mass flow rates were measured as part of the arc jet operations, while the measured 
pressures and heat fluxes at the stagnation point of ealorimeters (of various shapes and sizes) were used to determine 
the speeifie total enthalpy on the nozzle eenterline. More details of the measurements, ineluding details of 
ealorimetry, ean be found in the paper of Skokova.^ 

A complete set of experimental data is now available for use in simulation of the flow field that develops in the 
nozzle and over the test artiele. The simulation procedure was originally developed around v3.05 of the in-house 
flow solver, Dplr, and one of the objeetives of the present work is to upgrade the procedure with the latest version, 
v4.02.2,^^ which offers a more stable option of a subsonic inflow boundary condition based on the method of 
eharaeteristics. Detailed flow eomputations are performed first for the various ealorimeters used in the test. Onee 
the flow solver is ealibrated to calorimeter measurements (pressure and heat flux), detailed eomputations are 
performed for the test article. The resulting heat flux and pressure are then provided as inputs to v2.6.1 of the 
materials thermal response code, FiatP 

A. Flowfield Simulations 

The bulk enthalpy and total flow rate (along with mass fraetions of N 2 , O 2 , and Ar) of the test are provided as 
inputs to the previously referenced Nozzle _Throat_Conditions. This code, developed by Gokgen and Saunders,^^ 
uses the Cea eode^^ (formerly ealled the Gordon-McBride code) as the computational engine; the arc-heated gas 
mixture ean be assumed to be thermodynamie equilibrium on aceount of high pressure in the are eolumn. The 
resulting mixture mass density, velocity (obtained from the speed of sound and assumed Mach number), 
temperature, and mass fraetions of the eonstituent species (N 2 , O 2 , NO, N, O, and Ar) are used as pointwise inflow 
eonditions to Dplr. The computational domain consists of the entire convergent-divergent nozzle and the free jet in 
which a calorimeter or test artiele is inserted. The computational boundaries for the free jet (plus calorimeter) are 
truneated so that a supersonic extrapolation boundary condition can be applied to them. Further, the water-eooled 
nozzle wall and the surface of the copper calorimeter are assumed to be at a constant temperature of 400 K (127 °C), 
and fully catalytic to atom recombination. 

B. Simulations of Calorimeters: Round 1, Dplr v3.05 

In the first round of simulations, only calorimeters are considered. Computational grids for the ealorimeter 
geometries of AHF 295 had already been ereated and tested in pre-test computations. The grid block corresponding 
to eaeh of the ealorimeters was adapted to the bow shoek, and the wall-normal spacing controlled so as to aehieve a 
eell Reynolds number of 0(1) over the ealorimeter. The in-house code Sage^^ was used for the adaption of the grid 
bloek around the calorimeter. 

The computed values of pressure and cold-wall heat flux at the stagnation point of eaeh ealorimeter are given in 
Table 1. Also shown in the table are the effective radii of each non-hemispherical calorimeter, using the 4-inch 
hemispherical calorimeter as the reference. 
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Table 1. Computational predictions of pressure and cold-wall heat flux at the stagnation point using Dplr v3.05 




Experiment 



CFD {Dplr v 3.05), Z/buik 

= 13.7 MJ/kg 


Calorimeter 

/^stag 

^^stag 

Rifi 

/^stag 

^stag 

Diff. from expt. 

^eff‘ 


kPa 

W/cm^ 

in 

kPa 

W/cm^ 

Pressure 

kPa 

Heat flux 

W/cm^ 

in 

4in Hemi 

8.3 

222.0 

2.0 

7.99 

264.9 

-3.5% 

±19.3% 

2.0 

4in Iso-q 

8.4 

167.8 

3.501 

8.08 

190.1 

-3.6% 

±13.3% 

3.884 

4in Flat Face 

8.3 

123.0 

6.515 

8.14 

129.5 

-1.9% 

±5.3% 

8.369 

6in Flat Face^ 

8.5 

98.0 

10.263 

8.26 

102.7 

-2.8% 

±4.8% 

13.306 


^Arithmetic average of all measurements 
^Results for center slug only 

^Results using the hemispherical calorimeter as the reference in both experiment and CFD 


Although the agreement between eomputation and experiment is exeellent for the flat-faeed calorimeters, the 
disagreement for the calorimeters with surface curvature is large. This disagreement is even more disconcerting 
given that the experimental data appear to be internally consistent - they all yield nearly similar values for centerline 
enthalpy and effective radius computed for each non-hemispherical calorimeter (using the hemispherical one as the 
reference). 

Additional computations were performed for the hemispherical and iso-q calorimeters with flat-faced slugs 
(sensing surface) of 0.38-inch diameter, i.e., exposed faces of the copper slugs at the stagnation point were assumed 
to have zero curvature (or equivalently, an infinite radius of curvature). The pressure and heat flux distributions for 
the 4 calorimeters - hemispherical, iso-q, and flat-face cylindrical - are shown in Fig. 5. Also shown (as dashed 
lines) are the distributions for the hemispherical and iso-q calorimeter with flat-faced slugs. Clearly, the assumption 
of curvature mismatch, i.e., a flat- faced slug in a spherical section, does not help much to improve the accuracy of 
predictions. At best, the assumption of a flat- faced slug decreases the predicted heat flux by about 5%. The error 
bars (±5% for pressure, and ±15% for cold-wall heat flux) shown in the figure are somewhat arbitrary in that they do 
not represent actual measurement errors, and are based on estimates provided by experienced arc jet test engineers. 
They are only meant to provide a relative measure of calorimeter performance. 




Figure 5. Pressure and cold-wall heat flux distributions over various calorimeters employed in AHF 295. The 
symbols represent experimental data, but the ±15% error bars shown are assumed estimates of performance 
of a calorimeter, and meant to provide a basis for discussion. For pressure, the assumed error bars are ±5%. 


Based on these results, one is thus left with the uncomfortable situation of predictions in excellent agreement for 
some, but not all, calorimeters. Given the number of free parameters - inflow Mach number, wall temperature. 
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surface catalysis, etc. - an attempt has been made to upgrade the simulation procedure to the latest version (v4.02.2) 
oiDPLR, which has added new features. These results are discussed next. 

C. Simulations of Calorimeters: Round 2, Dplr v4.02.2 

Version 4.02.2 of Dplr offers a new feature of being able to specify a stagnation pressure and temperature for a 
subsonic inflow boundary, and uses the method of characteristics in conjunction with the assumption of isentropic 
flow. This should be contrasted with the procedure adopted in v3.05, where an inflow Mach number was specified 
and the inflow boundary was considered to be fixed (invariant in time). The inflow Mach number was determined 
approximately, and the Dirichlet boundary condition specifies more than the number of requisite variables, which 
could result in non-conservation of mass. 

As a first step towards handing off the process from v3.05 to v4.02.2 of Dplr, computations were performed for 
the 4-inch hemispherical calorimeter. The inflow conditions of v3.05 were used in v4.02.2, but the inflow Mach 
number was allowed to float. The resulting predictions of distributed surface temperature and cold-wall heat flux are 
shown in Fig. 6. 



(a) Pressure 



(b) Cold-wall heat flux 


Figure 6. Pressure and cold-wall heat flux distributions, obtained using v3.05 and v4.02.2 of Dplr for a 4-inch 
hemispherical calorimeter at AHF 295 conditions. 


The preliminary conclusion from this one-off computation is that the agreement in predictions of v3.05 and 
v4.02.2 of Dplr is fair. However, the reduction of one free parameter (inflow Mach number) in simulations is the 
primary motivator for wanting to switch to the newer version. In the simulation process based on v3.05 of Dplr, the 
inflow Mach number assumed at the nozzle inlet is usually 0.053 - based on simple application of gas dynamic 
equations, while the subsonic inflow boundary condition used in v4.02.2 yields a value of 0.0538, which is 1.5% 
higher. Despite the small difference in inflow Mach, one should realize that the Mach number is subsonic and has 
significant influence (~5%) on stagnation pressure (Fig. 6a). 

D. Simulations of Calorimeters: Round 3, Dplr v4.02.2, Catalycity effect study 

It should be noted that even with the elimination of the inflow Mach number as an uncertainty there is no 
improvement in the predicted level of heating at the stagnation point, and therefore another explanation should be 
sought. The surface temperature at the stagnation point of the test article (SPRITE-Tl) was monitored using IR 
pyrometers during Run 3 of AHF295. The time traces shown in Fig. 7 perhaps provide a clue. 
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(a) IR pyrometer located in test chamber (b) IR pyrometer located outside test chamber 

Figure 7. Time traces of surface temperature measured using (a) an IR pyrometer located inside the test 
chamber, and (b) a IR pyrometer located outside the test chamber during Run 3 of AHF 295. 

Time traces of test article surface temperature, acquired with an IR camera and a pyrometer, of the test article are 
between t ~ 1572 s to ^ ~ 1625 s (a period of roughly 50 s - the exposure time for the SPRITE test article), and are 
clearly seen in Fig. 7. The two smaller (and narrower) pulses on either side of the tall/broad pulse in Fig. 10 
represent the response of the pyrometers to the emission of the slugs of the 4-inch iso-q calorimeter inserted in the 
freestream before and after the exposure of the test article. The time traces are remarkably steady over the exposure 
times of the calorimeter, and both pyrometers record the elevated temperatures. 

It is hypothesized here that the copper slugs oxidize almost immediately upon insertion into the arc-heated 
stream, even if good care is taken in cleaning and polishing the slug prior to its use in arc jet measurements. The 
oxide layer that develops on the copper slug, perhaps several monolayers thick, is partially catalytic to 
recombination of atomic species diffusing to the surface of the calorimeter, i.e., the surface is not efficient enough to 
guarantee complete recombination of diffusing N and O atoms to their respective molecular states. It is hypothesized 
further that the oxide layer reaches radiative equilibrium very quickly and that its emissivity is similar to that of the 
test material (PICA). 

Based on the above hypothesis and from the pyrometer data, it is inferred that the exposed face of the slug 
reached about 800 °C (1073 K) during the few seconds that the calorimeter dwelled in the free jet, and that value 
was nearly a constant during the exposure. The value of approximately 840°C (1 1 13 K) read by the pyrometer inside 
the test chamber differs from that read by the pyrometer located outside, but that difference is not expected to impact 
heat computations performed with a temperature value of 800 °C (1073 K), which is larger than the nominal value 
of 127 °C (400 K) assumed in cold- wall computations. 

All results presented so far have assumed copper (material that makes up both the nozzle and the calorimeter) is 
fully catalytic to recombination of atoms. While this is a reasonable assumption to make if the wall is maintained at 
a sufficiently low temperature (such as the nozzle wall, which is water cooled), it becomes increasingly tenuous as 
the wall temperature exceeds 1000 K (the slug is not water cooled). A sufficiently convincing model for gas-surface 
interaction does not exist currently. However, the research of Barbato et al}^ can be put to practical use. As per the 
reference cited, the efficiencies of copper in recombining nitrogen and oxygen atoms are 0.28 and 0.1, respectively. 
For a slug surface temperature of 1073 K (and nozzle wall at 400 K), additional computations are performed for a 
noncatalytic slug surface (a limiting case), and a partially catalytic surface with recombination efficiencies of 
Barbato et al. The computational predictions of surface pressure and heat flux are shown in Fig. 8. 
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Figure 11. Influence of slug surface catalytic recombination efficiency on stagnation point pressure and heat 
flux for a 4-inch hemispherical copper slug calorimeter. 


The results shown in Fig. 8 indieate that eopper is likely to be partially eatalytie and definitely not non-eatalytie. 
While the predieted pressure is virtually insensitive to the surfaee eatalysis model employed, the heat flux for a 
partially eatalytie wall, or a wall with a redueed eatalytie effieieney to atom reeombination, is in mueh better 
agreement with experiment - predietion is only 4.2% higher than experiment. It then remains to be seen if the 
assumption of a partially eatalytie wall ean replieate measured pressure and heat flux for all ealorimeters used in 
AHF 295. The results of sueh eomputations are shown in Table 2. The laek of agreement in the ‘theoretieaf 
effeetive radius {i.e., effeetive radius inferred from CFD results) and that from the work of Zoby and Sullivan still 
remains to be investigated. It is likely that the boundary-layer edge eonditions are different for the different 
ealorimeters despite the faet that all of them were tested at the same distanee from the nozzle exit plane; flow 
relaxation between shoek and body depends on the shoek standoff distanee, whieh is largest for the 6-ineh flat-faeed 
eylindrieal ealorimeter and smallest for the hemispherieal one. 


Table 2. Computational predictions of pressure and cold-wall heat flux at the stagnation point using Dplr v4.02.2 


Calorimeter 

/^stag 

Experiment 

^stag 


/^stag 

CFD {Dplr v4.02.2), Ffb«ik = 13.7 MJ/kg 
^stag Diff. from expt. 



kPa 

W/cm^ 

in 

kPa 

W/cm^ 

Pressure 

kPa 

Heat flux 

W/cm^ 

in 

4in Hemi 

8.28 

222.0 

2.0 

8.34 

231.6 

+0.7% 

+4.2% 

2.0 

4in Iso-q 

8.38 

167.8 

3.501 

8.43 

169.1 

+0.6% 

+0.8% 

3.752 

4in Flat Face 

8.30 

123.0 

6.515 

8.49 

117.6 

+2.3% 

-4.4% 

7.757 

6in Flat Face 

8.50 

98.0 

10.263 

8.62 

94.1 

+1.4% 

-4.0% 

12.115 


^Results using the hemispherical calorimeter as the reference in both experiment and CFD 


The agreement between eomputation and experiment at the stagnation point is now exeellent for all ealorimeters. 
The pressure and heat flux distributions for the 4 ealorimeters - hemispherieal, iso-q, and flat-faee eylindrieal - are 
shown in Fig. 9. 
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Figure 9. Pressure and cold-wall heat flux distributions over various calorimeters employed in AHF 295. 
Computations assume copper to be partially catalytic to atom recombination, and the calorimeter surface is 
set to a constant temperature of 1073 K (pyrometer measurement at the stagnation point). 


The analysis presented thus far is a refinement over the CFD proeess built around v3.05 of Dplr. Although the 
analysis shows that it is quite plausible that copper is partially catalytic, further work and application to a larger 
number of test cases remains to be done. The good agreement between experiment and eomputations for this 
partieular are-heater setting for all ealorimeters is partieularly eneouraging. The lessons learned from simulations of 
ealorimeters are applied to the aetual test artiele. 

E. Test Article Simulation 

Simulation of the flow field around the test artiele is straightforward. The only additional ‘eomplexity’ in the 
eomputation is inelusion of the entire test artiele (sans the sting to avoid a time-eonsuming 3D eomputation), and 
part of the diffuser as well. The inflow eonditions are the same as those used in ealorimeter simulations with v4.02.2 
of Dplr. The assumption here is that the simulation proeedure has been suffieiently benehmarked against 
ealorimetry, whieh is foeused on the stagnation point alone. There are limited numbers of off- stagnation 
measurements with a 6-ineh flat-faeed ealorimeter. These data eover about 4 inehes (diameter) of the eore flow. 
Sinee the predieted heat flux at the off-stagnation loeations are in good agreement with measurements, making some 
allowanee for flow asymmetry, applieation of the eonditions and proeedures to the test artiele is likely going to be 
adequate. 

Computations have been performed using v4.02.2 of Dplr for the test artiele. As with the ealorimeters, 
eomputations have been performed three different ways: (a) the baseline method (designated as ewl) with nozzle, 
test artiele, and diffuser walls set to 400 K, and eatalytie reeombination effieieney set to unity, (b) the alternate 
method (designated as ew2) with the nozzle and diffuser walls set to 400 K, the test artiele surfaee set to 1073 K, 
and the eatalytie reeombination effieieneies (for N and O) set to values preseribed by the model of Barbato et al., 
and (e) the baseline hot wall method (designated as hw) in whieh the nozzle and diffuser wall are assumed fully 
eatalytie and set to a temperature of 400 K; the test article is assumed fully catalytic to atom recombination, but with 
the ability to re-radiate heat with an emissivity of 0.85. It should be noted that the third model does not eonsider 
eonduetion of heat through the thiekness of the material, i.e., an adiabatie baek wall is taeitly assumed in the model. 

The grid topology for these eomputations is not very amenable to grid tailoring or grid adaption, so no attempt is 
made here to adapt the grid to the bow shoek and boundary layer. Sinee the bow shoek will be misaligned with the 
grid, the heat flux, espeeially in the vieinity of the stagnation point, is expeeted to be a little less smooth. 
Furthermore, the eomputation does not inelude a sting, and the flow in the wake is likely to be unsteady 
(numerieally) past the flow separation point on the aft shell. 

In the absenee of a ealorimeter of the same shape and size as the test artiele, sealing arguments are relied upon to 
assess the aeeuraey of eomputations. Sinee the nose radius of the SPRITE-Tl geometry is 4-inehes, the eold-wall 
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heat flux should be roughly 41% lower than that for the hemisphere based on the Fay-Riddell eorrelation. From 
Table 2, this works out to roughly 155-160 W/em^ - the upper estimate coming from the fact that the nose radius of 
the SPRITE-Tl geometry and that of the 4-inch iso-q calorimeter are very similar (3.57 vs 4.00 inches). 

Distributions of surface pressure and heat flux obtained from the three computations mentioned above are shown 
in Fig. 10. Also shown are the predicted surface shear stresses and temperatures in Fig. 11. 




Figure 10. Pressure and heat flux (cold and hot wall) distributions over the SPRITE-Tl configuration tested 
in AHF 295. Results of computations in which all walls (nozzle, test article, and diffuser) are set to 400 K are 
designated by ‘cwl’ (blue curves) and computations in which the test article surface temperature is set to 
1073 K, and the nozzle and diffuser walls set to 400 K are designated by ‘cw2’ (green curves). The hot-wall 
computations are designated by ‘hw’ (red curves). 




(b) Surface temperature (adiabatic back wall) 


Figure 11. Shear stress and temperature distributions over the SPRITE-Tl configuration tested in AHF 295. 
Results of computations in which all walls (nozzle, test article, and diffuser) are set to 400 K (127 °C) are 
designated by ‘cwl’ (blue curves) and computations in which the test article surface temperature is set to 
1073 K (800 °C), and the nozzle and diffuser walls set to 400 K (127 °C) are designated by ‘cw2’ (green 
curves). The hot-wall computations are designated by ‘hw’ (red curves). 


The predicted radiation equilibrium temperatures (albeit with an adiabatic back wall) in the hot-wall 
computations are in excess of 2000 °C. At such temperatures, the assumption of a fully catalytic wall is tenuous. 
However, it is assumed that the heat transfer coefficient, required in materials thermal response computations, is 
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weakly dependent on wall temperature. Therefore, the results of the hot-wall eomputations for the SPRITE-Tl 
eonfiguration are used in response eomputations. 

F. Materials Thermal Response Analysis 

The eomputed flow field solution for the hot-wall case was post-proeessed using an in-house software utility, 
Blayer^^ to extraet boundary-layer properties - edge, wall, thiekness, ete. The surfaee pressure, edge enthalpy, and 
eonveetive heat transfer eoeffieient (film eoeffieient) were extraeted at loeations eorresponding to the eenter of eaeh 
instrumented plug (see Fig. 1). These quantities were provided as inputs to v2.6.1 of the materials thermal response 
eode. Fiat. The staek-up of the material on the heatshield was defined as PICA with 0.01 ineh thiek RTV-560 to 
bond the material to 0.125 thiek Al-2024 (the struetural shell). While the thiekness of PICA is 1 ineh over the 
eonieal frustum, it is about 1.36 inehes thiek at the stagnation point, beeause the aluminum shell had a flat faee 
rather than a eurved one eonforming to the nose radius of the OML. Eaeh instrumented plug had 3 in-depth K-type 
thermoeouples. The 2 plugs on the eonieal frustum had the thermoeouples at nominal depths of 0.125, 0.375, and 
0.625 ineh, and for the stagnation point plug these depths were 0.25, 0.50, and 0.75 ineh. The inputs used in Fiat 
eomputations are shown in Table 3. Computations were run out to 3600 seeonds, whieh ineluded 50 seeonds of 
exposure of the test artiele to the are-heated stream. 


Table 3. Inputs to Fiat\2,6.\ 


Plug location 

PICA 

thickness 

in 

TTT fiber 
orientation 

H 

BTU/lbm 

Pe^e^H 

Ibm/tf .s 

P 

atm 

Stagnation point 

1.36 

0° 

5869 

0.0310 

0.0828 

Mid-cone 

1.00 

45° 

5877 

0.0126 

0.0379 

Bottom cone 

1.00 

45° 

5878 

0.0104 

0.0377 


Predieted time histories of surfaee temperature and thermoeouple temperatures are shown in Fig. 12 through Fig. 
14, respeetively for the three instrumented plugs at the apex, mid-frustum, and end- frustum regions of the test 
artiele. In these figures. Fiat predietions are shown as red lines, and measurements (pyrometer or thermoeouple 
data) are shown as blaek lines. 




Time, s Time, s 

(a) Surface temperature (b) In-depth temperature 

Figure 12. Time histories of (a) Fiat predicted surface temperatures (red lines) and measured surface 
temperatures (black lines), and (b) Fiat predicted temperatures in the material at 3 depths from the surface 
and thermocouple measurements. Results are shown for the instrumented plug at the apex (spherical nose 
cap) of the model (see Fig. 1). 
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Time, s Time, s 

(a) Surface temperature (b) In-depth temperature 

Figure 13. Time histories of (a) Fiat predicted surface temperatures (red lines) and measured surface 
temperatures (black lines), and (b) Fiat predicted temperatures in the material at 3 depths from the surface 
and thermocouple measurements. Results are shown for the instrumented plug at the middle of the conical 
frustum (see Fig. 1). 




(a) Surface temperature (b) In-depth temperature 

Figure 14. Time histories of: (a) Fiat predicted surface temperatures (red lines) and measured surface 
temperatures (black lines), and (b) Fiat predicted temperatures in the material at 3 depths from the surface 
and thermocouple measurements. Results are shown for the instrumented plug at the end of the conical 
frustum (see Fig. 1). 


The surface temperature trends appear to match well against predictions (Figs. 15a, 16a, and 17a), as do the 
trends at the in-depth thermocouple locations. At a time instant of / = 25 s (mid-point of the 50 s exposure time for 
the test article), values of temperature measured using pyrometers are compared with predicted surface temperatures 
in Table 4. Predictions are always above the measurements by as much as 80 °C. Whether these differences are due 
to higher inflow enthalpy has not yet been explored fully. Further, for the instrumented plug at the center of the 
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conical frustum, there is a spike in temperature at the time instant that the model is taken out of the are-heated 
stream. It is speculated that cracks that appeared in the vieinity of the instrumented plugs could be the problem.^ 


Table 4. Comparison of measured and predicted surface temperatures 


Plug location 

Experiment (@t=25 s) 

Computation (@t=25 s) 

AJ=Je.p-7’c«mp 


°C 

°C 

°C 

Stagnation point 

1962 

2044 

-78 

Mid-cone 

1513 

1559 

-46 

Bottom cone 

1405 

1468 

-63 


For the in-depth thermocouples, the level of agreement between measurements and predictions depends on the 
measurement location. For instance, at the stagnation region (the apex of the test artiele) the predieted rise in 
temperature (Fig. 12b) traeks that of eaeh thermocouple. However, predictions indicate more rapid eool down past 
the 50 s mark (the end of exposure of the model to the arc-heated stream). For the plug at the middle of the conical 
frustum, the agreement is mueh better than at the stagnation point, although the model still prediets more rapid 
eooling (Fig. 13b) than is observed. For the instrumented plug towards the shoulder of the test artiele, the predieted 
temperatures still show rapid deereases (Fig. 14b) eompared to measurements, after the model is taken out of the 
are-heated stream. Fiat is a one-dimensional material response eode, and whether two-dimensional effects are 
important has not been explored here. However, the paper of Skokova et alJ explores the use of a 2D eharring 
material ablation program called Titan It is also speeulated that the proximity of the model to the are-heated 
stream, despite being moved out of it, eould have an effeet. 

Version 2.6.1 of permits analysis of sensitivity of materials response to variations in enthalpy. In the present 
work a ±10% variation (around the nominal) in enthalpy is considered. The predieted recession for the various cases 
is shown in Table 5. The results indicate a linear dependence of reeession on enthalpy - a ±10% ehange in enthalpy 
yields a ±12% change in reeession. 


Table 5. Predicted recession using Fiat\2,6.\ 



Stagnation Point Plug 

Mid-cone Plug 

Bottom-cone Plug 


Recession 

in 

Difference 

Recession 

in 

Difference 

Recession 

in 

Difference 

0.9 Nom. enthalpy 

0.1855 

-10.8% 

0.0658 

-11.5% 

0.0523 

-11.9% 

Nom. enthalpy 

0.2079 


0.0743 


0.0594 


1.1 Nom. enthalpy 

0.2305 

10.8% 

0.0829 

11.5% 

0.0664 

11.9% 


More details of reeession measurements ean be found in the work of Skokova et alJ 

V. Concluding Remarks 

The CFD-based process for analyzing the flow of an arc-heated mixture expanded through a eonvergent- 
divergent nozzle was used to determine the maximum diameter of 45° sphere-cone models (DS2 shape or SPRITE- 
T1 shape) that could be tested safely without flow bloekage in the 18-inch nozzle of the AHF at NASA ARC. Based 
on the eomputations, a 14-ineh diameter wooden model of the DS2 was built and tested suceessfully in the AHF. 
Apart from proving that there would be no issues with flow bloekage for a model of this size, the charred region of 
the back portion of the wood model provided a visual confirmation of flow separation. The length of the region of 
attaehed flow (indicated by eharring) on the baek shell eorrelated well with CFD simulations. 

The eomputational proeess was then applied to the 14-ineh diameter SPRITE-Tl model assuming nominal 
eonditions (elose to the maximum operational eapability of the AHF), and the predieted aerothermal environments 
were used to perform thermal- struetural analysis. The thermal analysis provided an estimate of the time of exposure 
of the test artiele so that the temperature of the battery of the data aequisition system built into the probe would not 
exceed 50 °C. 

The computational process was used in simulating the flow around calorimeters of various sizes and shapes of 
test AHF 295. Actual arc-heater settings were used in the study of the flow around ealorimeters. One uneertainty 
that has not been resolved in the present work is the anomalous values^ of bulk enthalpy measurements made using 
the enthalpy by energy balanee or EB2 method.^^ The average value of 11.5 MJ/kg reported by the EB2 
measurements is significantly lower than data gathered over several years in the AHF, and mueh lower than the 13.7 
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MJ/kg estimate for the enthalpy on the nozzle eenterline. This differenee suggests that the inflow to the eonvergent- 
divergent nozzle is non-uniform, a finding that is not borne out by a eombination of heat flux-pitot sweeps of the 
free jet. This diserepaney in EB2 measurements was ignored in the present work, and eomputations performed 
assuming a uniform inflow with a bulk enthalpy of 13.7 MJ/kg. 

The eomputational proeess, developed around v3.05 of Dplr, was upgraded to one based on v4.02.2. This 
upgrade utilizes the new eapability of Dplr to handle subsonie inflows with a eharaeteristies based method, thus 
eliminating the need to speeify an inflow Maeh number (as in the older framework). However, it should be noted 
that this gain is eurrently limited to uniform inflow eonditions, and Dplr will require further modifieations to handle 
pointwise speeifieation of inflow mass flow and enthalpy. The resulting predietions were found not to have as mueh 
internal eonsisteney as the measurements, and it was hypothesized that the eopper slugs used in the ealorimeters 
were partially catalytic to atom (N and O) reeombination at the surfaee. Using a model of eopper eatalyeity 
developed by Barbato et al. flow simulations were performed again, and this time the predieted results were not only 
in good agreement with measurements aeross all shape/size ealorimeters used in the test. Despite its sueeess, the 
hypothesis of partial eatalyeity remains to be tested aeross several operating eonditions of both the AHF and IHF 
before any definitive eonelusions ean be made. 

Next, the CFD-based proeess was applied to the aetual test artiele, SPRITE-Tl, assuming the entire model 
(ineluding the aft shell) to be fully eatalytie to atom reeombination. The predieted aerothermal environments 
(notable pressure and heat flux) were used in the material thermal response eode. Fiat. The results of Fiat 
predietions were eompared against surfaee temperature measurements made with optieal pyrometers, and with 
temperature traees of in-depth thermoeouples plaeed in the forebody TPS material, PICA. The eomparisons were 
found to be fair for surfaee temperatures. For in-depth temperatures, while the temperature rise was well predieted 
for all thermoeouples for the exposure period of 50 s. Fiat predieted a more rapid eool down than measurements. 
Whether this has something to do with material eraeking (observed after taking the model out of the are-heated 
stream), or with two-dimensional heat eonduetions effeets, or even proximity to the heated stream after the model 
has been taken out, still remains to be investigated. Having a ealorimeter of the same size and shape as the SPRITE- 
Tl test artiele would have also helped in determining whether there were flow non-uniformities in the test. 
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